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ABSTRACT  Magnesium  is  crucial  for  baculovirus  transmission  in  Culex  nigripalpus  (Theobald)  and 
Culex  quinquefasciatus  (Say)  larvae.  However,  the  mechanistic  role  of  magnesium  in  baculovirus 
transmission  is  unknown.  To  investigate  the  possible  role  of  host  response  factors  in  baculovirus 
transmission,  suppression  subtractive  hybridization  was  used  to  identify  genes  differentially  tran¬ 
scribed  after  magnesium  exposure  in  Cx.  quinquefasciatus  larvae.  Suppression  subtractive  hybridiza¬ 
tion  was  performed  in  both  directions  enriching  for  cDNAs  differentially  transcribed  between  a 
nonmagnesium  larval  control  and  magnesium  ( 15  mM  MgS04)  treatment  of  Cx.  quinquefasciatus  larvae 
held  for  1  h  at  27°C.  Clones  from  differentially  transcribed  genes  were  evaluated  by  sequencing,  and 
relative  gene  transcription  levels  were  analyzed  using  quantitative  real-time  polymerase  chain  reac¬ 
tion  quantitative  real-time  polymerase  chain  reaction.  Target  transcripts  up/downregulated  by  mag¬ 
nesium  included  Cx.  quinquefasciatus  troponin  C,  isocitrate  dehydrogenase,  allergen,  cytochrome  b5, 
chymotrypsinogen,  apolipophorins,  tryptase  gamma,  carboxylesterase,  prolylcarboxypeptidase,  imag- 
inal  disc  growth  factor,  aldehyde  dehydrogenase,  tropomyosin-1,  chitotriosidase,  heat  shock  protein 
70  B2,  inorganic  phosphate  cotransporter,  and  many  other  hypothetical  protein  genes.  Magnesium  can 
alter  gene  transcription  in  a  vector  mosquito  population,  and  understanding  this  process  can  provide 
insight  into  the  mechanistic  role  of  magnesium  in  baculovirus  transmission, 
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The  critical  role  of  divalent  ions  in  the  larval  medium 
for  the  mediation  (Mg^^)  and  inhibition  (Ga^^,Gu^^, 
Fe^^,  and  Zn^^)  of  Culex  nigripalpus  (Theobald)  del- 
tabaculovirus  (GuniNPV)  infection  in  Cufer  larvae  has 
been  reported  (Becnel  et  al,  2001).  Infections  of  Gun¬ 
iNPV  are  restricted  to  nuclei  of  midgut  epithelium, 
with  mortality  occurring  72-96  h  postinfection.  Mag¬ 
nesium  is  only  required  for  the  initial  stage  of  the 
infection  process  before  the  virions  entering  the  mid- 
gut  cells  (Becnel  et  al.  2001).  The  high  infectivity  and 
pathogenicity  of  GuniNPV  against  the  principal  vec¬ 
tors  of  West  Nile  virus  in  North  America  make  Gun¬ 
iNPV  an  attractive  candidate  for  future  development 
as  a  biopesticide  (Becnel  2006),  Metal-induced  gene 
expression  in  the  mosquito  midgut  has  been  previ¬ 
ously  documented:  exposing  Aedes  aegypti  (L.)  larvae 
to  heavy  metals  showed  AEIMUGl  (Ae.  aegypti  intes- 
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tinal  mucin)  to  be  a  metal-responsive  gene,  MRG 
(Rayms-Keller  et  al.  2000),  and  molecular  biological 
studies  have  identified  midgut  genes  induced  by  heavy 
metals  in  the  midgut  epithelium  of  aquatic  arthropods 
(Beaty  et  al,  2002). 

In  the  current  study,  suppression  subtractive  hy¬ 
bridization  (SSH)  is  used  to  identify  genes  that  are 
differentially  transcribed  during  magnesium  treat¬ 
ment  of  larval  Culex  quinquefasciatus  (Say).  Expres¬ 
sion  patterns  are  confirmed  by  quantitative  real-time 
polymerase  chain  reaction  (qPGR) ,  and  the  regulated 
genes  that  were  identified  are  discussed.  In  addition, 
a  time  course  study  of  gene  transcription  in  response 
to  Mg^^  in  Culex  larvae  was  also  examined  for  known 
genes.  Characterization  of  genes  induced  by  Mg^^  in 
Culex  larvae  could  provide  a  mechanistic  understanding 
of  gene  expression  in  larvae  exposed  to  heavy  metals. 

Materials  and  Methods 

RNA  Extraetion.  Cx.  quinquefasciatus  ( Orlando,  FL 
strain,  maintained  since  1952)  were  reared  in  the  in¬ 
sectary  of  the  Mosquito  and  Fly  Research  Unit  at  the 
Center  for  Medical,  Agricultural,  and  Veterinary  En¬ 
tomology,  United  States  Department  of  Agriculture- 
Agricultural  Research  Service  (Gainesville,  FL)  fol¬ 
lowing  previously  established  protocols  (Gerberg  et 
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al.  1994).  Magnesium  treatment  (15  mM  MgS04)  of 
third  instar  Cx.  quinquefasciatus  larvae  (250-300)  was 
carried  out  for  1  h  at  27°C  for  the  polymerase  chain 
reaction  (PCR) -select  cDNA  subtraction.  This  mag¬ 
nesium  concentration  has  been  shown  to  be  optimal 
for  CuniNPV  transmission  (Becnel  et  al.  2001).  Con¬ 
trol  Cx.  quinquefasciatus  larvae  (250-300)  were  col¬ 
lected  at  room  temperature  (27°C).  The  total  RNAs 
were  extracted  using  TRIzol  reagent,  according  to  the 
manufacturer’s  instructions  (Invitrogen,  Carlsbad, 
CA).  Poly  (A)  ^  RNA  was  isolated  applying  Oligo- 
tex-dT  suspension  (Qiagen,  Valencia,  CA).  RNA  and 
mRNA  samples  were  quantified  by  SmartSpec  Plus 
spectrophotometry  (Bio-Rad,  Hercules,  CA). 

Subtractive  Hybridization.  Differentially  expressed 
mRNA  populations,  both  magnesium  (Mg^^)  treat¬ 
ment  and  control,  were  converted  into  cDNA  using 
the  PCR-Select  cDNA  Subtraction  Kit  (BD  Clontech, 
Roche  Molecular  Systems,  Alameda,  CA) .  Subtractive 
cloning  was  performed,  as  previously  described,  and 
differentially  expressed  cDNA  sequences  were  used 
to  construct  a  subtracted  cDNA  library  (Patel  and  Sive 
2001,  Zhao  et  al,  2009) .  In  the  generalized  subtraction 
scheme,  the  mosquitoes  to  be  compared  are  the  [  -I-  ] 
or  tracer  (Mg^^  treatment)  and  the  [  — ]  or  driver 
(non-Mg^^  treatment,  control),  in  which  niRNAs  ex¬ 
pressed  in  the  tracer  and  not  the  driver  are  isolated. 
Hybrids  that  form  include  sequences  common  to  both 
mosquito  populations.  The  unhybridized  fraction  is 
enriched  for  sequences  that  are  preferentially  ex¬ 
pressed  in  the  tracer  mosquito  population.  Differen¬ 
tially  expressed  cDNA  sequences  were  used  to  con¬ 
struct  a  subtracted  cDNA  library  for  further  study. 

PCR-Select  cDNA  Subtraction  Library.  Forward 
and  reversed  subtracted  libraries  were  cloned  using 
the  TOPO  TA  Cloning  Kit  for  sequencing  (Invitrogen, 
Carlsbad,  CA).  Transformed  plasmids  were  inserted 
into  One  Shot  TOPIC  Competent  Cells  (Invitrogen) 
and  grown  overnight  on  Luria-Bertani  (LB)  plates 
containing  ampicillin  and  X-gal  (5-bromo-4-chloro-3- 
indolyl-jS-D-galactopyranoside).  When  the  cell  trans¬ 
formation  was  conducted,  the  cells  were  incubated  in 
the  37°C  shaker  for  exactly  I  h,  and  plated  out  in  the 
LB  +  ampcillin  -I-  X-gal  plates.  In  that  case,  we  limited 
the  clone  number  in  the  plate  (slightly  over  100)  so  as 
they  were  not  overly  crowded.  We  picked  out  almost 
all  clones  on  the  plates.  For  each  library,  >100  white 
colonies  were  isolated  and  grown  overnight  in  LB- 
ampicillin  broth  at  37°C  and  235  rpm  in  the  Innova 
4000  Incubator  Shaker  (New  Brunswick  Scientific, 
Edison,  NJ). 

Gene  Sequencing  of  PCR-Select  cDNA  Subtraction 
Library.  Clones  of  the  subtracted  library  were  purified 
with  QIAprep  Minipreps  (Qiagen,  Valencia,  CA) .  The 
0.5  pg  of  plasmid  DNAs  was  then  digested  by  using  2.5 
U  of  EcoRI  enzyme  for  1.5  h,  and  was  run  on  a  1% 
agarose  gel  to  confirm  the  DNA  insert.  Selected  clones 
were  then  sequenced  by  the  Interdisciplinary  Center 
for  Biotechnology  Research,  University  of  Florida 
(Gainesville,  FL).  The  ABI  Big  Dye  terminator  se¬ 
quencing  was  used  for  DNA  sequencing.  The  se¬ 
quences  were  analyzed  using  the  National  Center  for 


Biotechnology  Information  BLASTN  program  to  iden¬ 
tify  sequence  homologies.  The  homologous  DNA  frag¬ 
ments  in  the  cDNA  or  mRNA  were  also  recorded  in 
Table  1.  Hypothetical  and  conserved  hypothetical 
proteins  were  further  analyzed  using  the  ExPASy  Web 
site  (http:/ /ca.expasy.org)  “DAS”  Transmembrane 
Prediction  Server  to  predict  potential  transmembrane 
segments. 

cDNA  Synthesis  for  qPCR.  A  3-p.g  aliquot  of  puri¬ 
fied  RNA  (from  «“I00  larvae  per  sample)  was  reverse 
transcribed  in  a  20-p,l  reaction  volume  using  Clone 
AMV  First-Strand  Synthesis  Kit  and  using  oligo(dT)2o 
primer  for  cDNA  synthesis,  according  to  the  manu¬ 
facturer’s  instructions  (Invitrogen,  Carlsbad,  CA).  To 
ensure  that  no  genomic  DNA  contaminated  the  sam¬ 
ple,  we  used  only  oligo(dT)  primer  to  synthesize  the 
first-strand  cDNA  library  for  the  reverse  transcription 
PCR.  The  reaction  was  terminated  by  heat  inactiva¬ 
tion  at  95°C  for  5  min.  The  cDNA  samples  for  qPCR 
from  Mg^^  treatment  and  control  were  diluted  by 
adding  80  pi  of  ddHaO  (300  ng/pl)  and  stored  at 
-20°C. 

To  design  gene-specific  primers,  detailed  analyses 
of  the  nucleotide  sequence  of  genes  found  in  the 
library  were  performed  using  the  PRIMER3-Design 
Primer  Pairs  and  Probes  program  from  Biology  Work¬ 
bench  (http:/  / workbench.sdsc.edu) .  The  primers  for 
the  Cx.  quinquefasciatus  actin  gene  (GenBank  acces¬ 
sion  number:  XM  001866828)  were  also  designed  for 
internal  control  and  comparison.  The  primers  are 
listed  in  Table  2, 

qPGR  Amplification.  The  qPCR  assay  for  Culex 
genes  was  performed  using  Platinum  SYBR  Green 
qPCR  SuperMix-UDG  with  BOX  (Invitrogen,  Carls¬ 
bad,  CA)  in  a  volume  of  15  pi  on  an  Applied  Biosys¬ 
tems  7300  Fast  Real-Time  PCR  System  (Foster  City, 
CA),  The  PCR  mixture  consisted  of  I  pi  of  diluted 
cDNA  (300  ng/pl),  0,5  pM  primers,  and  IX  master 
mix.  In  every  qPCR  run,  actin  was  used  as  an  internal 
control  to  normalize  for  variation  in  the  amount  of 
cDNA  template.  The  PCR  primers  for  actin  gene  were 
CxActin629-F  5'-TGCGTGACATCAAGGAGAAG-3' 
and  CxActin-893R  5'-GTGTTGGCGTACAGGTC- 
CTT-3'.  The  PCR  thermal  cycling  parameters  were 
the  same  as  previously  described  (Zhao  et  al.  2008, 
2009).  Relative  expression  levels  were  calculated  as 
follows.  First,  CxGENE  transcript  levels  relative  to  a 
standard  (CxACTIN)  were  calculated  using  the  for¬ 
mula  A  cycle  threshold  (CT)  =  CT  (CxGENE)  —  CT 
(CxACTIN).  Second,  an  average  ACT  value  for  each 
sample  was  calculated.  Third,  relative  expression  lev¬ 
els  were  calculated  using  the  equation  100  X  2 
[  —  average  ACT  ] .  Lastly,  fold  increases  or  decreases 
were  calculated  using  Mg^^  treatment  RNA  relative 
expression  level/non-Mg^^  treatment  control  RNA 
relative  expression  level.  The  detailed  analyses  of  the 
nucleotide  sequence  of  all  up/downregulated  genes 
were  compared  between  Mg^^  treatment  and  control. 
When  the  relative  gene  expression  was  >1,  the  gene 
expression  of  treatment  was  upregulated.  However, 
when  the  relative  gene  expression  was  <1,  the  gene 
expression  of  treatment  was  downregulated. 
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Table  1.  Genes  identified  from  Culex  quinquefasciatus  larvae  in  response  to  magnesimn  exposure 


Clone 


Accession  no.,  organism,  putative  identity,  mPlNA/cDNA  length 


Known  proteins 
lAB-3,  lAB-16, 
lAB-53 
lAB-5 

lAB-14.  lAB-35, 
lB-9 
lAB-18 
lAB-ia  lB-6 
lAB-25 
lAB-36 
lAB-42 
lAB-51 
lB-4 

lB-8,  lB-17 
lB-10 

lB-20 

lB-24 

lB-29 

5B-7 

Hypothetical  proteins 
lAB-4 
lAB-15 
lAB-29 
lAB-34 
lAB-44 
lB-25 

Conserved 

hypothetical 
proteins 
lB-7,  lB-14 
lB-21 
lB-22 
lB-32 
5B-5 

Ribosomal  proteins 
lAB-8,  lAB-9, 
lB-53,  lAB-52 


refXM_001844392.1 

refXM_001844978.1 

refXM_001869433.1 

refXM_001863979.1 
refXM_001861639.1 
refXM_001849258.1 
refXM_001847441.1 
refXM_001868370.1 
ref  XM  001850434.1 
refXM_001844621.1 
ref  XM  001845150.1 
refXM_001841858.1 
length  =  1,407 
refXM_001864200.1 
refXM_001867666.1 
ref  XM  001861401.1 
ref  XM  001847074.1 
length  =  1,569 

refXM_001851302.1 
refXM_001858971.1 
refXM_001862348,l 
refXM_001847371.1 
refXM_001850204.1 
refXM  001847210.1 


refXM_001866472.1 

refXM_001851203.1 

refXM_001859612,l 

refXM_001848166.1 

refXM_001861615.1 

refXM  001864716.1 


Cx.  quinquefasciatus  troponin  C  partial  mRNA,  length  =  1,421 

Cx.  quinquefasciatus  isocitrate  dehydrogenase  partial  mRNA,  length  =  1,303 
Cx.  quinquefasciatus  allergen  partial  mRNA,  length  =  911 

Cx.  quinquefasciatus  cytochrome  b5  partial  mRNA,  length  =  554 
Cx.  quinquefasciatus  chymotrypsinogen  partial  mRNA,  length  =  801 
Cx.  quinquefasciatus  apolipophorins  partial  mRNA,  length  =  10,169 
C.r.  quinquefasciatus  tryptase  gamma  pai*tial  mRNA,  length  =  828 
Cx.  quinquefasciatus  carboxylesterase  partial  mRNA,  length  =  1,695 
C.t.  quinquefasciatus  prolylcarboxypeptidase  partial  mRNA,  length  =  1,506 
C.r.  quinquefasciatus  imaginal  disc  growth  factor,  mRNA,  length  =  1,384 
Cx.  quinquefasciatus  aldehyde  dehydrogenase,  mRNA,  length  =  1,499 
C.r.  quinquefasciatus  isocitrate  dehydrogenase  cytoplasmic,  mRNA, 

Cr.  quinquefasciatus  tropomyosin- 1,  mRNA,  length  =  875 
C.r.  quinquefasciatus  chitotriosidase-1,  mRNA,  length  =  2,643 
Cr.  quinquefasciatus  heat  shock  protein  70  B2,  mRNA,  length  =  1,917 
Cr.  quinquefasciatus  inorganic  phosphate  cotransporter,  mRNA, 


Cr.  quinquefasciatus  hypothetical  protein  partial  mRNA,  length  =  1,068 
Cr.  quinquefasciatus  hypothetical  protein  partial  mRNA,  length  =  1,311 
Cr.  quinquefasciatus  hypothetical  protein  partial  mRNA,  length  =  827 
Cr.  quinquefasciatus  hypothetical  protein  partial  mRNA,  length  =  558 
Cr.  quinquefasciatus  hypothetical  protein  partial  mRNA,  length  =  2,078 
Cr.  quinquefasciatus  hypothetical  protein,  mRNA,  length  =  2,000 


Cr.  quinquefasciatus  conserved  hypothetical  protein,  mRNA,  length  =  747 
Cr.  quinquefasciatus  conserved  hypothetical  protein,  mRNA,  length  =  657 
Cr.  quinquefasciatus  conserved  hypothetical  protein,  mRNA,  length  =  612 
Cr.  quinquefasciatus  conserved  hypothetical  protein,  mRNA,  length  =  1,998 
Cr.  quinquefasciatus  conserved  hypothetical  protein,  mRNA,  length  =  691 

Cr.  quinquefasciatus  ribosomal  protein  S25  partial  mRNA,  length  =  566 


lAB-11,  lAB-40 

lAB-21 

lAB-24 

1AB-.32,  lB-12 

lAB-49 

lB-13,  lB-28 

lB-23 

lB-27 


gbDQ401 444.1  Cr.  quinquefasciatus  clone  CpiLD-AOl  large  subunit  ribosomal  RNA,  partial 
sequence;  mitochondrial,  length  =  231 

ref  XM_001861287.1  Cr.  quinquefasciatus  60S  ribosomal  protein  L14  partial  mRNA,  length  =  660 
ref  XM_001848606.1  Cr.  quinquefasciatus  60S  ribosomal  protein  L3  partial  mRNA,  length  =  1,423 
ref  XM_001846289.1  Cr.  quinquefasciatus  60S  ribosomal  protein  L29  partial  mRNA,  length  =  455 
ref  XM_001842837.1  Cr.  quinquefasciatus  60S  ribosomal  protein  L32  partial  mRNA,  length  =  613 
ref  XM_001841875.1  Cr.  quinquefasciatus  60S  ribosomal  protein  L8,  mRNA,  length  =  998 
refXM_001845506.1  Cr.  quinquefasciatus  60S  ribosomal  protein  L4,  mRNA,  length  =  1,477 
ref XM_001863381.1  Cr.  quinquefasciatus  ribosomal  protein  Lll,  mRNA,  length  =  815 


Gene  regions 


nt;  1146-767,  nt: 

1187-720 
nt:  504-1057 
nt:  911-395,  nt: 

239-878 
nt:  281-113 
nt:  801-535 
nt:  10122-9807 
nt:  75-381 
nt:  321-493 
nt:  1328-1506 
nt:  435-733 
nt:  927-1274 
nt:  951-1290 

nt:  875-827 
nt:  2484-2217 
nt:  821-1287 
nt:  177-1 


nt:  1040-495 
nt:  300-1 
nt:  639-127 
nt:  270-498 
nt:  735-858 
nt:  1269-1711 


nt:  510-708 
nt:  1-102 
nt:  612-292 
nt:  1830-1993 
nt:  685-501 

nt:  159-490,  nt: 

490-30,  nt: 

32-490 
nt:  231-61 

nt:  373-602 
nt:  945-393 
nt:  18-349 
nt:  218-27 
nt:  636-801 
nt:  579-274 
nt;  25-428 


Time  Course  Study.  Magnesium  treatment  (15  mM 
MgS04)  of  third  instar  Cx.  quinquefasciatus  larvae  was 
carried  out  for  15  min,  1  h,  2  h,  and  4  h  at  27°C  («=250 
larvae  each  time  point)  to  trace  the  RNA  expression 
level  of  differentially  transcribed  genes  as  identified 
from  SSH.  Control  Cx.  quinquefasciatus  larvae  (^250 
larvae  each  time  point)  were  simultaneously  collected 
without  magnesium  treatment  at  room  temperature 
(27°C).  The  time  course  study  was  replicated  three 
times.  The  RNA  was  extracted  and  cDNA  was  created, 
as  previously  described  in  this  study,  and  qPCR  was 
performed  using  primers  designed  from  genes  found 
by  SSH  (Tables  1  and  2). 

Statistical  Analysis.  Comparisons  of  means  were  an¬ 
alyzed  using  Student’s  t  test,  and  t  values  and  P  values 
were  reported  when  normality  and  equal  variance 
tests  were  passed.  Significant  differences  between  the 


data  were  determined  using  SigniaPlot  software 
(SigmaPlotll.2,  Systat  Software,  San  Jose,  CA). 

Results 

Identification  of  Genes  Specifically  Transcribed 
During  Magnesium  Treatment  in  Cx.  quinquefascia¬ 
tus.  To  reveal  how  larval  genes  were  specifically 
transcribed  during  magnesium  exposure  in  Cx.  quin¬ 
quefasciatus,  we  used  subtractive  hybridization  to 
construct  a  subtraction  cDNA  library  from  control 
larvae  and  larvae  that  had  been  exposed  to  magne¬ 
sium  for  f  h.  Clones  of  differentially  transcribed 
genes  were  examined  by  sequencing.  Using  National 
Center  for  Riotechnology  Information  blast  to  an¬ 
alyze  the  sequences,  we  found  22  clones  encoding 
Cx.  quinquefasciatus  known  genes  representing  15 
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Table  2.  Primers  used  for  quantitative  RT-PCR 


Primer  name“  sequence 

1AB-3R  5'-ATGATGACCGGlTAAGCTGC-3' 
lAB  .3F  5'-AAAGGGGAGGGAGTITCTTG-3' 
1AB-4R  5'-GAGACCGAAGAAAAACACGC-3' 
1AB-4F  5'-TTGGTTGGTTGGGAGAATGT-3' 
lAB-llR  5'-CGACGTCGATGTTGGArrAAG-3' 
lAB-llF  5'-AGAAACGAAGGTGGCTrAGG-3' 
1AB-15R  5'-GGTCCTGAAACACTCGGAGA-3' 
1AB-15F  5'-GAAGnTCGAGGTCCCTTCC-3' 
1AB-18R  5'-GTGAGGAGGTGGTGGTTGA-3' 
1AB-18F  5'-AGGGCGCTTTGTrCATmT-3' 
1AB-21R  5'-CGGTGGCGTTAACGGTTAAC-3' 
1AB-21F  5'-GAAGAAGGAGGGGAAGAAGA-3' 
lAB-24R5'-ATGGACCATGTGGTCCTGTC-3' 
lAB-24F5'-GAGATCGGTGTGATGAAGCA-3' 
1AB-25R  5'-GGGGGAAGGGATAAACTGTA-3' 
1AB-25F  5'-TAATGGGGTGTTGGATGATG-3' 
1AB-29R  5'-GAGTGAGTGAGCGAGGAAGC-3' 
1AB-29F  5'-TATCCGCTCACAATTCCACA-3' 
1AB-34R  5'-GGACATCTTCTTGAAGGCCA-3' 
lAB  .34F  5'-GTACGGCAAGAACCGGAAG-3' 
1AB-36R  5'-TGAATGCGATCGATGGTAAA-3' 
1AB-36F  5'-AGCCTGGGTTATCTCGTTGA-3' 
1AB-42R  5'-GCTGAAGAACCCAAGCAGAC-3' 
1AB-42F  5'-GGTGCGGAATTAAGGATGGT-3' 
1AB-44R  5'-CTTCGTCGACGTCAGACTTG-3' 
1AB-44F  5'-AGACGAAGCCCCGAATGTAT-3' 
1AB-49R  5'-ATGGAAnTGTAGCCAAGGG-3' 
1AB-49F  5'-AAGAAGGGCAAGCTGATTGT-3' 
1AB-51R  5'-TGATGAATGTGGTGGTrTGC-3' 
1AB-51F  5'-GATCGAGGGCAATGGTrATT-3' 
1B-4R  5'-AGTCCATCGACGTTGnTCC-3' 

1B-4F  5'-CCCGCCATCATGAACTATCT-3' 

1B-6R  5'-ACG'mTCGTrCAGAlTGGG-3' 

1B-6F  5'-GTTACTGCGGGTCATGTGCT-3' 

1B-9R  5'-GTTCTGCGACATGGTGTGAT-3' 

1B-9F  5'-ACCTAAGCTGAAACGGGGAT-3' 

IB-IOR  5'-CAGACCTTCTCCAGCGTr-3' 

IB-lOF  5'-AGTACCAGCAGGGCAAAGAA-3' 
1B-12R  5'-GAGAAAAGGGCATGTTGAAA-3' 

1B-12F  5'-AATGTGCCAAAAGTrCGTGG-3' 

1B-13R  5'-ACCGGTCmTCCTCCAGAT-3' 

1B-13F  5'-ACTTCCGTGACCCGTACAAG-3' 
1B-14R  5'-GTGGTGCAACTTAGGGGACT-3' 
1B-14F  5'-AGCTGAGCTGAAGAAGGTGG-3' 
1B-17R  5'-TGTrGTCCGTAACATGCGAA-3' 

1B-17F  5'-GGATGAGGAGCAGTTCAACA-3' 
1B-20R  5'-TAGCGACGAGGTTGAACTCC-3' 

1B-20F  5'-CAGTGCGGTGATGAGAACAA-3' 
1B-21R  5'-GGTACTTGGCATGGTGGTGT-3' 

1B-21F  5'-TGGACGCTAGTATGGGGAAG-3' 
1B-22R  5'-CTGGGGTGTGGGTATGTTGT-3' 

1B-22F  5'-TGGCGTTAACTCGACGTTCT-3' 

1B-23R  5'-CGTGGTrCTGACGGTTGTTC-3' 

1B-23F  5'-ATGAAATCGGGAAGGTTCTG-3' 

1B-24R  5'-GTGGTGGTCGTTGTAGGG-3' 

1B-24F  5'-GTGGGGAAATCGTAGCAAGT-3' 
lB-25R5'-ATGGTGGTGGTGTGCCTAGT-3' 

1B-25F  5' TrAAAACCGTGACACGGTGA-3' 

1B-27R  5'-AAGCAGTrCATrGGGTCCTG-3' 
lB-27F5'-GGTCGATGGGTATGTAGGGT-3' 

1B-29R  5'-ATCAACTrCGTGAAGGGTGG-3' 

1B-29F  5'-GGGAGAAGGACGAGTTrGAG-3' 
1B-32R  5'-GTGTGCTTGGTrGGAATGGT-3' 

1B-32F  5'-GAGGATATrGTGCTGGTTGG-3' 
lB-53R5'-TGTCGTAGGTGGGCTTATCG-3' 

1B-53F  5'-CGGAAAAGACGCAGAAGAAG-3' 

5B-5R  5  -ATCGTTGGCGTAGTrGCTCT-3' 

5B-5F  5'-GTGAACCGCTGGTGAAAATG-3' 

5B-7R  5'-CGACGACAGAATGTACCCCT-3' 

5B-7F  5'-AAGAGTATCTGGCGGAGGGT-3' 

GxActm-629F5'-TGCGTGACATGAAGGAGAAG-3' 

CxActin-893R5'-GTGTrGGCGTACAGGTCCTT-3' 


“  Primer  names  correspond  to  the  clone  numbers  in  Table  1. 


different  genes:  troponinC,  isocitrate  dehydroge¬ 
nase,  allergen,  cytochrome  b5,  chymotrypsinogen, 
apolipophorins,  tryptase  gamma,  carboxylesterase, 
prolycarboxypeptidase,  imaginal  disc  growth  factor, 
aldehyde  dehydrogenase,  tropomyosin-1,  chitotrio- 
sidase-1,  heat  shock  protein  70B2,  and  inorganic 
phosphate  cotransporter.  We  also  found  11  clones 
encoding  six  hypothetical  and  five  conserved  hypo¬ 
thetical  proteins.  In  addition,  we  found  15  clones  that 
encoded  for  Cx.  quinquefasciatus  cytosolic  ribosonial 
subunit,  60S  ribosomal  proteins,  and  one  mitochon¬ 
drial  large  subunit  ribosomal  RNA  (Table  1). 

Magnesium  Effects  on  RNA  Relative  Expression 
Levels  of  Different  Genes  in  Larval  Cx.  quinquefas¬ 
ciatus.  qPCR  analyses  were  conducted  to  examine  the 
differential  transcription  of  genes  identified  in  re¬ 
sponse  to  magnesium  in  Cx.  quinquefasciatus  larvae. 
According  to  our  qPCR  data,  most  genes  present  in  the 
subtracted  cDNA  library  were  upregulated  (Fig.  lA). 
Compared  with  the  control,  the  RNA  relative  expres¬ 
sion  level  of  apolipophorins  (lAB-25),  prolylcar- 
boxypeptidase  (lAB-51),  imaginal  disc  growth  factor 
(lB-4),  and  tryptase  (lAB-36)  increased  significantly 
10-fold  for  1-h  Mg^^  treatment  of  Cx.  quinquefasciatus 
larvae  (Fig.  lA;  Supplemental  Table  3A  [available 
online  only]).  Troponin  (lAB-3,  lAB-16,  and  lAB- 
53),  allergen  (lAB-14,  lAB-35,  lB-9),  cytochrome  b5 
(llAB-18),  chymotrypsinogen  (lAB-19,  and  lB-6), 
carboxylesterase  (lAB-42),  aldehyde  dehydrogenase 
(lB-8,  and  lB-17),  isocitrate  dehydrogenase  (lB-10), 
tropomyosin-1  (lB-20),  chitotriosidase-1  (lB-24),and 
heat  shock  protein  70  B2  (lB-29)  increased  «“2-  to 
6-fold  for  the  1-h  Mg^^  treatment  of  Cx.  quinquefas¬ 
ciatus  larvae  (Fig.  lA).  Inorganic  phosphate  cotrans¬ 
porter  (5B-7)  found  in  the  reverse  subtractive  cDNA 
library  was  downregulated  for  the  1-h  Mg^^  treatment 
of  Cx.  quinquefasciatus  larvae  (Fig.  lA).  All  of  these 
genes  identified  in  response  to  Mg^"*"  exposure  for  1  h 
were  significantly  different,  when  compared  with  the 
control  (Supplemental  Table  3  [available  online 
only]).  A  >9-fold  increase  in  the  transcription  of 
genes  of  hypothetical  proteins  (e.g.,  lAB-4,  lAB-15, 
and  lAB-44)  and  conserved  hypothetical  protein 
inRNA  (lB-14  and  lB-21)  was  found  after  the  1-h 
Mg^^  treatment  of  Cx.  quinquefasciatus,  which  was 
significantly  different  from  the  control  (Fig.  IB,  Table 
1;  Supplemental  Table  3B  [available  online  only] ). 

In  addition,  RNA  expression  levels  of  several  ribo¬ 
somal  genes  were  modulated  and  upregulated  more  than 
two  times.  The  qPCR  analysis  further  substantiated  the 
differential  transcription  of  the  mitochondrial  large  sub¬ 
unit  ribosomal  gene  (GenBank  accession  number; 
gb  DQ401444.1)  in  response  to  Mg^^  exposure  (Fig.  IC; 
Supplemental  Table  3C  [available  online  only] ). 

We  also  conducted  a  time  course  study  of  magne¬ 
sium  exposure  on  the  relative  RNA  expression  levels 
of  several  different  genes  in  larval  Cx.  quinquefasciatus 
to  better  understand  early  gene  response  and  regula¬ 
tion.  RNA  expression  levels  of  many  genes  were  mod¬ 
ulated,  some  showing  a  >5-fold  increase  in  expression 
during  the  Mg^^  exposure.  For  example,  the  imaginal 
disc  growth  factor  gene  was  significantly  upregulated 
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lAB-3  lB-9  lAB-18  IB-B  lAB-25  lAB-3e  lAB-42  lAB-51  lB-4  lB-17  lB-10  lB-20  lB-24  lB-29  5B-7 


lB-53  lAB-21  lAB-24  lB-12  lAB-49  lB-13  lB-23  lB-27  lAB-11 

Fig.  1.  The  qPCR  showed  the  relative  gene  transcription  (fold  increase)  after  magnesium  treatment  for  1  h  compared 
with  the  control.  (A)  Known  genes  (Table  1)  differentially  transcribed  after  magnesium  treatment  for  1  h  compared  with 
the  control  (nonmagnesium  treatment).  (B)  Unknown /Conserved  genes  (Table  1)  differentially  transcribed  after  magne¬ 
sium  treatment  for  1  h  compared  with  the  control.  (C)  Ribosomal  genes  (Table  1)  differentially  transcribed  after  magnesium 
treatment  for  1  h  compared  with  the  control. 


>10-fold  after  1  h  of  magnesium  exposure  (Fig.  2A; 
Supplemental  Table  4A  and  G  [available  online 
only] ) .  This  same  gene,  however,  was  downregulated 
after  4  h  of  magnesium  exposure  (Fig.  2A) .  According 
to  our  qPCR  data,  relative  gene  transcription  of  car- 
boxylesterase,  chitotriosidase,  and  apolipophorin  have 
similar  patterns  (i.e.,  upregulated  at  both  15  min  and 
1  h  of  Mg^^  exposure,  and  downregulated  after  2  and 
4  h  of  Mg^^  exposure)  (Fig.  2B,  C,  and  F;  Supple¬ 
mental  Table  4B,  C,  F,  and  G  [ available  online  only] ) . 
For  chymotrypsinogen,  relative  gene  transcription 
was  downregulated  after  4  h  of  Mg^^  exposure  (Fig. 


2D;  Supplemental  Table  4D  and  G  ]  available  online 
only]).  Although  the  relative  gene  transcription  of 
cytochrome  b5  was  upregulated  throughout  the  time 
course  of  Mg^^  exposure,  it  reached  *=8.7-fold  more 
than  the  control  after  1-h  Mg^^  exposure  and  was 
significantly  different  (Fig.  2E;  Supplemental  Table  4E 
and  G  ] available  online  only] ). 

Discussion 

Genes  upregulated  in  response  to  elevated  magne¬ 
sium  concentrations  could  be  involved  in  salt  regula- 
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Fig.  2.  The  qPCR  showed  the  relative  gene  transcription  (fold  increase)  during  the  time  course  (15  min,  1  h,  2  h,  and 
4h),  (A)  Relative  imaginal  disc  growth  factor  gene  (lR-4)  transcription  (fold  increase)  after  magnesium  treatment  compared 
with  the  control  at  15  min,  1  h,  2  h,  and  4  h.  (R)  Relative  carboxylesterase  gene  (lAR-42)  transcription  (fold  increase)  after 
magnesium  treatment  compared  with  the  control  at  15  min,  1  h,  2  h,  and  4  h.  (C)  Relative  chitotriosidase  gene  (lR-24) 
transcription  (fold  increase)  after  magnesium  treatment  compared  with  the  control  at  15  min,  1  h,  2  h,  and  4  h,  (D)  Relative 
chymotrypsinogen  gene  (lR-6)  transcription  (fold  increase)  after  magnesium  treatment  compared  with  the  control  at  15  min, 
1  h,  2  h,  and  4  h.  (E)  Relative  cytochrome  b5  gene  (lAR-18)  transcription  (fold  increase)  after  magnesium  treatment 
compared  with  the  control  at  15  min,  1  h,  2  h,  and  4  h.  (F)  Relative  apolipophorin  gene  (lAR-25)  transcription  (fold  increase) 
after  magnesium  treatment  compared  with  the  control  at  15  min,  1  h,  2  h,  and  4  h. 


lion  and  other  physiological  functions  in  mosquitoes, 
but  might  also  serve  to  enhance  baculovirus  infection 
in  the  midgut.  It  has  been  reported  that  magnesium  is 
crucial  for  transmission  of  baculovirus  CuniNPV  in  Cx. 
nigripalpus  and  Cx.  quinquefasciatus  (Becnel  et  al. 
2001).  It  was  also  shown  that  calcium  is  a  potent 
inhibitor  of  transmission  of  CuniNPV  even  in  the  pres¬ 
ence  of  magnesium.  Our  current  study  shows  that 
magnesium  levels  that  are  sufficient  for  CuniNPV 
transmission  can  alter  the  gene  transcription  of  larval 


Cx.  quinquefasciatus.  There  are  15  known  genes  and  1 1 
unknown  genes  that  are  up/  downregulated  during  the 
Mg^^  exposure  in  Culex  larvae  (Table  1,  Figs.  1  and  2) , 
Specifically,  genes  that  are  overexpressed  in  the  mid¬ 
gut,  where  baculovirus  entry  and  replication  occur, 
may  play  a  role  in  enhanced  viral  transmission. 

Carboxylesterases  are  ubiquitous  enzymes  and  play 
diverse  metabolic  roles  in  mammalian  organisms,  in¬ 
sects,  and  plants,  and  in  the  detoxification  of  a  broad 
range  of  drugs  and  xenobiotics  (Fonnum  et  al.  1985; 
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Potter  and  Wadkins  2006;  Cui  et  al.  2007a,  2007b; 
Cummins  et  al.  2007;  Holmes  et  al.  2009).  One  study 
showed  that  the  activities  of  carboxylesterases  in  re¬ 
sponse  to  three  insecticides  (i.e.,  malathion,  propoxur, 
and  permethrin)  differed  by  mosquito  species;  high 
activities  of  carboxylesterases  were  detected  in  Culex 
tritaeniorhynchus  (Giles),  but  not  Culex  gelidus 
(Theobald)  (Karunaratne  and  Hemingway  2000) .  In  a 
previous  study,  the  expression  of  carboxylesterase 
genes  in  the  mosquito  Culex  pipiens  (L.)  was  impor¬ 
tant  in  resistance  to  organophosphate  insecticides 
(Cui  et  al,  2007a).  In  Ae.  aegypti  larvae,  induction  of 
several  cytochrome  P450  genes  and  a  carboxylesterase 
gene  by  exposure  to  low  doses  of  various  xenobiotics 
(e.g.,  insecticide,  copper)  increased  larval  tolerance 
to  insecticides  (Poupardin  et  al.  2008).  These  authors 
suggested  that  exposure  of  mosquito  larvae  in  natural 
sites  to  low  levels  of  xenobiotics  affects  tolerance  to 
insecticides,  possibly  through  cross-induction  of  cer¬ 
tain  detoxification  genes.  In  this  study,  there  was  a 
significant  3.2-fold  increase  in  the  transcription  of  a 
carboxylesterase  gene  induced  by  1-h  magnesium  ex¬ 
posure.  Induction  of  a  carboxylesterase  gene  by  mag¬ 
nesium  may  also  increase  the  tolerance  of  larval  mos¬ 
quitoes  to  pathogens  by  some  unknown  mechanism 
and  could  have  an  impact  on  virus  infection  in  the 
midgut  of  Culex  larvae. 

Recently,  mammalian  chitinases  have  attracted 
much  attention,  because  chitotriosidase  plays  a  key 
role  in  the  pathogenesis  of  asthma  and  in  plasma  of 
Gaucher  patients  (Hollak  et  al.  1994;  Boot  et  al.  1995; 
Renkema  et  al.  1995, 1998;  Zhu  et  al.  2004).  Chitotrio¬ 
sidase  secreted  by  activated  human  macrophages  is 
involved  in  the  defense  against  chitin-containing 
pathogens  such  as  fungi,  nematodes,  and  insects 
(Grosso  et  al.  2004,  Gordon-Thomson  et  al.  2009). 
Studies  by  Di  Luca  et  al.  (2006,  2007)  showed  that 
human  chitotriosidases  from  different  sources  helped 
Plasmodium  falciparum  (Welch)  enter  the  Anop/iefes 
midgut  presumably  through  degradation  of  the  peri- 
trophic  matrix  (PM).  According  to  our  data,  chito¬ 
triosidase  gene  transcription  was  slightly  upregulated 
after  15-min  (1.47-fold)  or  1-h  (4.13-fold)  M^^  ex¬ 
posure  compared  with  the  control  and  downregulated 
after  2-  and  4-h  Mg^^  exposure  (Fig.  2C;  Supplemental 
Table  4C  and  G  [available  online  only]).  Virions  of 
CuniNPV  must  pass  through  the  larval  PM  to  attach  to 
and  enter  midgut  cells,  and  this  occurs  during  the  first 
hours  of  exposure.  Degradation  of  PM  in  larval  mos¬ 
quitoes  by  chitotriosidase  produced  in  response  to 
magnesium  exposure  could  serve  to  enhance  infec- 
tivity  by  CuniNPV  in  the  midgut  of  Culex  larvae. 

Chymotrypsinogen,  a  precursor  of  the  digestive  en¬ 
zyme  chymotrypsin,  is  converted  to  chymotrypsin  by 
trypsin.  Chymotrypsinogen  is  involved  in  the  solubi¬ 
lization  and  activation  of  insecticidal  protein  toxins 
produced  by  Bacillus  thuringiensis  (Zhu  et  al.  1997) .  In 
this  study,  the  transcription  of  a  chymotrypsinogen 
gene  was  upregulated  (5.8-fold)  after  15  min  (Fig. 
2D),  The  virions  of  CuniNPV  are  encapsulated  by  an 
occlusion  body  protein  with  4-8  virions  per  occlusion 
body  (OB).  Once  OBs  are  ingested  by  mosquito  lar¬ 


vae,  the  occlusion  body  protein  must  be  solubilized  to 
release  virions  before  crossing  the  PM  and  attachment 
to  microvilli  of  midgut  cells  (Becnel  et  al.  2001 ) .  A  host 
response  to  Mg^^  exposure  that  increases  the  produc¬ 
tion  of  digestive  enzymes  such  as  chymotrypsinogen 
could  mediate  the  release  of  CuniNPV  virions  from 
OBs  and  plays  a  role  in  enhancing  virus  infection. 

Apolipophorins  contain  the  major  component  of 
lipophorin,  which  mediates  transport  for  various  types 
of  lipids  in  hemolymph  and  acts  by  forming  lipopro¬ 
tein  particles  that  bind  lipoproteins  and  lipids.  There 
are  two  mature  apolipophorins:  apolipophorin-I  and 
apolipophorin-II  (Marinotti  et  al.  2006).  In  the  insect 
hemolymph,  the  complex  of  low-density  lipophorin, 
high-density  lipophorin,  and  apolipophorin  III  pro¬ 
vides  a  reusable  lipid  shuttle  for  flight  muscle  energy 
supply  (Surholt  et  al.  1992).  An  apolipophorins  gene 
was  expressed  in  fat  body  tissues  throughout  devel¬ 
opment  in  Anopheles  gambiae  (Marinotti  et  al.  2006), 
In  a  previous  study,  the  expression  of  apolipophorins 
in  Ae.  aegypti  was  induced  by  blood  feeding,  and 
lipophorin  was  synthesized  in  the  fat  body  (van  Heus- 
den  et  al.  1998) .  Our  data  show  that  an  apolipophorins 
gene  increased  >10-fold  after  1-h  Mg^^  exposure 
(Fig.  2F).  There  were  also  significant  differences  in 
the  transcription  of  the  apolipophorins  gene  between 
magnesium  exposure  at  1  h  and  the  control  in  Culex 
larvae  (Supplemental  Table  4F  and  G  [available  on¬ 
line  only]).  The  transcription  of  apolipophorins  in 
larval  Cx.  quinquefasciatus  was  induced  by  Mg^^  ex¬ 
posure,  which  could  elevate  lipophorin  synthesis  in 
the  fat  body  and  change  the  lipid  structure  of  the  cell 
membrane. 

Some  hypothetical  proteins  were  highly  upregu¬ 
lated  in  response  to  Mg^^  exposure  (lAB-15,  lAB-4, 
lB-14,  lB-21,  and  lAB-44),  According  to  ExPASy 
(“DAS,”  Transmembrane  Prediction  Server),  all  of 
these  genes  have  potential  transmembrane  segments, 
but  their  possible  role  in  salt  regulation  and/ or  bac- 
ulovirus  transmission  is  unknown.  The  imaginal  disc 
growth  factor  (lB-4)  is  also  highly  upregulated,  but 
how  this  protein  functions  in  mosquitoes  is  unknown. 
Notably,  some  MRGs  found  in  previous  studies  with 
Diptera  such  as  AEIMUCI  (Ae.  aegypti  intestinal  mu¬ 
cin)  and  CTTUBI  {Chiromomus  tentans  a-tubulin) 
were  not  identified  in  this  study  (Beaty  et  al.  2002). 
Sublethal  concentrations  of  toxic  metals,  such  as  cop¬ 
per  and  calmium,  induced  expression  of  these  two 
proteins,  which  could  enhance  survivorship  of  aquatic 
arthropods  in  polluted  environments.  Because  mag¬ 
nesium  is  not  normally  toxic  to  mosquito  larvae, 
the  lack  of  overexpression  of  homologues  is  not 
unexpected. 

The  current  study  suggests  that  certain  genes  (e.g., 
chymotrypsinogen,  chitotriosidase,  carboxylesterases, 
apolipophorin)  expressed  in  response  to  magnesium 
exposure  may  play  a  role  in  the  interaction  of  Cx. 
quinquefasciatus  larvae  and  infection  by  the  baculo- 
viruses  CuniNPV,  Identification  of  these  genes  is  the 
first  step  in  determining  the  mechanistic  role  of  mag¬ 
nesium  in  transmission  of  baculovirus  in  mosquitoes, 
and  also  provides  basic  information  on  genes  that  may 


1060 


Journal  of  Medical  Entomology 


Vol.  47,  no.  6 


be  involved  in  other  physiological  responses  in 
mosquitoes. 
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